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Abstract 

The sum of the CP- violating asymmetries A(A°_) and in hyperon nonleptonic decays is presently 

being measured by the E871 experiment. We evaluate contributions to the asymmetries induced by 
chromomagnetic-penguin operators, whose coefficients can be enhanced in certain models of new physics. 
Incorporating recent information on the strong phases in 3 — > Air decay, we show that new-physics con- 
tributions to the two asymmetries can be comparable. We explore how the upcoming results of E871 may 
constrain the coefficients of the operators. We find that its preliminary measurement is already better than 
the e parameter of K-K mixing in bounding the parity-conserving contributions. 
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I. INTRODUCTION 



The phenomenon of CP violation remains one of the least understood aspects of particle physics. 
Although CP violation has now been detected in several processes in the kaon and P-meson sys- 
tems [1], its origin is still far from clear. The standard model (SM) can accommodate the experi- 
mental results, but they do not yet provide critical tests for it [2]. To pin down the sources of CP 
violation within or beyond the SM, it is essential to observe it in many other processes. 

Nonleptonic hyperon decays provide an environment where it is possible to obtain additional 
observations of CP violation. Although this has been recognized for a long time [3], only recently 
has it been experimentally feasible to search for CP violation in some of these decays [4-6]. There 
is currently such an effort being done at Fermilab, by the HyperCP (E871) Collaboration [5, 6]. 

The reactions studied by HyperCP are the decay sequence H~ — > An~ , A — > pn~ and its 
antiparticle counterpart. For each of these decays, the decay distribution in the rest frame of the 
parent hyperon with known polarization w has the form 

dr 

+ (1) 

where p is the unit vector of the daughter-baryon momentum and a is the parameter relevant to the 
CP violation of interest. By evaluating the decay chain S — > An — > pnn, the HyperCP experiment 
is sensitive to the sum of CP violation in the S decay and CP violation in the A decay. Thus it 
measures [5, 6] 

A SA = A A + A s , (2) 

where _ _ 

A* = A(EZ) = as + ® E , A A ee A(A°) ee ^A±^a (3) 

a s — a s a A - a A 

are the CP- violating asymmetries in S — > An and A — * pn, respectively. The published mea- 
surements currently available [7] are A A = 0.012 ± 0.021 and A SA = 0.012 ± 0.014. HyperCP will 
obtain more precise results, with an expected sensitivity of ~10~ 4 , and has recently reported [6] a 
preliminary measurement of A SA = (—7 ± 12 ± 6.2) x 10 -4 . 

The amplitudes for both A — > pn~ and 2~ — > An~ contain S- and P-wave components, each 
of which consists of contributions describing |AJ| = | and | transitions, with the former being 
known to dominate. Assuming |A7| = \ dominance, one can derive at leading order [8] 

A A = — tan (6p - 6s) sin (<j)p - <f>g) , A s = - tan (6$ - 5|) sin (0| - 0|) . (4) 

Here, 6g (6p) is the strong S- wave (P-wave) Nir scattering phase-shift at = M A , and 6^ (6p) is 
the strong 5-wave (P-wave) An scattering phase-shift at ^fs = M s . Moreover, (0p = ) are the 
CP-violating weak phases induced by the \AS\ = 1, |A/| = | interaction in the S*-wave (P-wave) 
of the A — > pn~ and 2~ — > A7r _ decays, respectively. 

The strong Nn scattering phases needed in Eq. (4) have been measured [9] to be 5$ ~ 6° and 
6p ~ —1° with errors of about 1°. In contrast, the strong An phases are less well determined. Using 
the current PDG numbers [7], one can deduce [10] the phase difference 6p — 6cj = (—7.7 ± 7.7)°. 
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Very recently, the E756 Collaboration [11] has published a new measurement of Sp — Sg = (+3.17± 
5.28 ± 0.73)°. HyperCP is presently also measuring this quantity, with better precision, and has 
reported [10] a preliminary result of 5p — 5$ = (7.6 ± 1.3^2 s) • 

On the theoretical side, the most recent update [12] of the standard-model prediction of A EA 
yields a value that is smaller than most of earlier estimates [8, 13], but with a sizable uncertainty, 
resulting in the range \A EA \ <5x 10~ 5 , which is compatible with some of the earlier predictions. 
Thus, the upcoming data from HyperCP will likely be insensitive to SM effects. Beyond the SM, 
the asymmetry is potentially more detectable, as various estimates [14, 15] indicate that A A could 
be as large as 10~ 3 in models with enhanced chromomagnetic-penguin operators (CMO). In these 
studies, the corresponding value of A E has been neglected because most of recent calculations based 
on chiral perturbation theory [16, 17] suggest that 5f> — 5| is small compared to 5p~ — 5$. However, 
there are early indications, from a coupled-channel K- matrix estimate given in Ref. [17] and from 
the preliminary result of HyperCP above [10], that the two phase-differences may actually be 
comparable in size. This is also compatible with the other two measurements of Sf, — 5| mentioned 
earlier. Since HyperCP is sensitive only to the sum A A + A E , it is therefore important to have an 
up-to-date expectation of A A s from possible new physics and of their sum. 

In this paper, we estimate both A A and A E due to possible physics beyond the SM, incorporating 
the new information on the strong phases and taking into account constraints from kaon data. 
Specifically, we consider contributions generated by the CMO, which could be significantly larger 
that their SM counterparts [14, 15]. The relevant effective Hamiltonian can be written as [18] 

H w = C g Q g + C g Q g + H.c. , (5) 
where C g and C g are the Wilson coefficients, and 

Q 9 = ij^ d ^ (1 + 7 5 ) * <%, , Q g = ^ doT (1 - 7b ) s G; v , (6) 

are the CMO, with being the gluon field-strength tensor and Tr(t a t b ) = \5 ah . These operators 
also contribute to the CP- violating parameters e in kaon mixing and e' in kaon decay, as well as 
to other hyperon and kaon observables [18-20]. Although e, e', and A A S receive contributions from 
the same \AS\ = 1 interaction, they probe different parts of it. Whereas e and e' are sensitive only 
to parity-even and parity-odd contributions, respectively, A A 3 are sensitive to both. Thus, with e 
and e' now being well measured, we will estimate the range of A EA arising from the CMO that is 
allowed by e and e', and then compare it with the preliminary result of HyperCP mentioned above. 
Since various new-physics scenarios may contribute differently to the coefficients of the operators, 
we will not focus on specific models, but will instead adopt a model-independent approach, only 
assuming that the contributions are potentially sizable. Accordingly, we will also explore how well 
the coefficients can be constrained in the event that HyperCP detects no CP-violation. 

In Sec. II, we employ heavy-baryon chiral perturbation theory to derive the decay amplitudes at 
leading order. In Sec. Ill, we calculate the weak phases using matrix elements estimated in the MIT 
bag model. We then estimate the CP- violating asymmetries, taking into account constraints from 
CP violation in the kaon system, and present a discussion of our results. We give our conclusions 
in Sec. IV. 
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II. CHIRAL LAGRANGIAN AND DECAY AMPLITUDES 



The chiral Lagrangian that describes the interactions of the lowest-lying mesons and baryons is 
written down in terms of 3 x 3 matrices tp and B which contain the lightest meson-octet and baryon- 
octet fields, respectively [21, 22]. The mesons enter through the exponential E = £ 2 = exp(i<^//), 
where / is the pion-decay constant. 

In the heavy-baryon formalism [22], the baryons in the chiral Lagrangian are described by 
velocity-dependent fields, B v . For the strong interactions, the chiral Lagrangian to lowest order 
in the derivative expansion is given by [22] 

CP = (B v iv ■ VB V ) + 2D (B V S» {A^ B V }) + 2F (B V S» [A^ B v ] > + \f (d^ d^) , (7) 

where (• • • ) denotes Tr(- • • ) in flavor-SU(3) space, S v is the spin operator, VB V = d»B v + [V, B v ] , 
with = § (f + g d^) , and A^ = \ (f - £ f d^) . In this Lagrangian, D and F are 
free parameters which can be determined from hyperon semileptonic decays. We will adopt the 
parameter values obtained from fitting tree-level formulas [22], namely D = 0.80 and F = 0.50. 
We will also need the chiral Lagrangian that explicitly break chiral symmetry [23], containing one 
power of the quark-mass matrix M = diag(m u , m d , m s ) , 

4 2) = \f (X + ) + ^ (B v { X+ ,B V }) + ^- (B v [ X+ ,B V }) + A- ( X+ ) (B V B V ) , (8) 

where we have used the notation x+ — £*x£*+£x*£ to introduce coupling to external (pseudo) scalar 
sources, x = s + ip, such that in the absence of the external sources x reduces to the mass 
matrix, x = 2B M. We will take the isospin limit, m u = m d = rh, and consequently x = 
diag(m 2 , m 2 ^ , 2m 2 K — m 2 ). In Eq. (8), the constants B , b DF0 are free parameters which can be 
fixed from data. 

In the weak sector, the chiral Lagrangian induced by the chromomagnetic-penguin operators 
has to respect their symmetry properties. From Eq. (5), we observe that Q g and Q g transform 
as (3 L ,3 R ) and (3 L ,3 R ), respectively, under SU(3) L xSU(3) R transformations. Accordingly, the 
desired Lagrangian at leading order is 

Av = P D (B v {?h?, B V }) + F (B v [?h?, B v ] > + P (htf) (B V B V ) + ^ f 2 B {htf) 

+ P D (B v {th£, B v }) + f3p (B v [£ht, B v }) + & (hT) (B V B V ) + ^ f 2 B (hE) + H.c. , (9) 

where (5 i are parameters containing the coefficient C g (C g ), and the 3x3-matrix h selects out 
s — > d transitions, having elements h kl = S k2 S 3l . As shown in Appendix A, the expression for C w 
can be inferred from the lowest-order chiral realization of the densities d(l ±7 5 )s. We remark that 
this Lagrangian is of 0(1) in the derivative and m s expansions. 

For the weak decay of a spin-| baryon B into another spin-^ baryon B' and a pseudoscalar 
meson 0, the amplitude in the heavy-baryon approach has the general form [24] 

iM B _ BI(t> = -i< J B'0|£ w+s | J B> = u B , (Af BI4> + 2S V - P4> A { B ^) u B , (10) 
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where the superscripts refer to the S- and P-wave components of the amplitude. These components 
are related to the decay width T and parameter a by 

r = 7 -\ e b> +m B ,) { \s \+ \p \) , a = | | 2 , | | 2 , ( n ) 

4:7im B ' \sy + \p\ z 

where s = A^ and p= \p B ,\A (p) . To express our results, we also adopt the notation [24] 



a 



(s,p) = k f A (S,P) 



ins', - v'2 / -4 . (12) 

From the Lagrangians given above, one can derive the S- and P-wave amplitudes at leading 
order, represented by the diagrams in Figs. 1 and 2, respectively. For the S'-wave, the first diagram 
is directly obtained from a weak vertex provided by Eq. (9), whereas the other diagram involves 
a ^-vacuum tadpole from Eq. (9) and a strong B^B'Kn vertex, which consists of contributions 
from both C^p and d 2 \ It is worth mentioning here that the and /3 terms in C w do not 
contribute to K — > mr decay, as the corresponding direct and tadpole diagrams cancel exactly [25]. 
For A — > pn~ and S~ — > , the resulting amplitudes are then 1 

(13) 



J s ) - _L (Q- _ 3 0-) _ JL 0- m z- m A 



m 



where (5 i = f3 i — (5 i and we have used the relations 

m A -m N = ~l(b D + 3b F ){m s - m) , m E - m A = l(b D - 3b F ){m s - m) , 

= P (m s + m) , 

derived from Eq. (8). For the P-wave, the amplitude arises from two baryon-pole diagrams, each 
involving a weak vertex from Eq. (9) and a strong vertex from Eq. (7), and a kaon-pole diagram 
involving a strong vertex from Eq. (7) followed by a K-ir vertex from Eq. (9). Thus we find 

fl(P) = (D + P) (gg + 3gg + 2P (gg - gg + QD + 3P) g+ 
Ap7r \/6 (m A — %) \/6 (m s — %) \/6 (m s — m) 

(15) 



(P) (-P + P) (gg - 3gg 2P (gg + /?+) (P-3P)/5+ 

a S-A7T- ~ ^ ' ^ ^ ' ^ 



(m s — m A ) \/6 (m H — m s ) \/6 (m s — m) 

where f3f = f3 i + f3 i and we have used m K — m 2 = B (m s — rh). We note that the baryon and 
meson masses in all the amplitudes above are isospin-averaged ones. 



1 The contribution from d 1 ^ to the tadpole amplitude contains the factor v ■ = m B — m B , = 0(m s ). As a result, 
the d 1 ^ and contributions to the amplitude both have the same m s order, C(m°), as the j3~ terms indicate. 
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FIG. 1: Leading-order diagrams for chromomagnetic-penguin contributions to S- wave hyperon nonleptonic 
decays. In all figures, a solid (dashed) line denotes a baryon (meson) field, and a solid dot (hollow square) 
represents a strong (weak) vertex. 



FIG. 2: Leading-order diagrams for chromomagnetic-penguin contributions to P-wave hyperon nonleptonic 
decays. 

III. NUMERICAL RESULTS AND DISCUSSION 



In order to estimate the weak phases in A A s , we first need to determine the parameters and P i 
in terms of the underlying coefficients C g and C g , respectively. From the effective Hamiltonian in 
Eq. (5) and the chiral Lagrangian in Eq. (9), we can derive the one-particle matrix elements 



/ \nj i A \ Pp + 3/? F + ~P D + 3/3 F _ 
{n\n w \A) = -j= u n u A 



<A|^ W |S°> 



Pd-ZPf + Pd-Wf - 

Vb A Uz 



(16) 



\/6 

+ p v ) B . 



Since there is presently no reliable way to calculate these matrix elements from first principles, we 
will employ the MIT bag model to estimate them, following earlier work [15]. Using the results 
given in Appendix B, and setting C = 0, we find 



P. 



D 



TlR 2 9 



-8V2I M N* 



M ' s "' l<_ Q 



ttB R 2 



where N, R, and I M are bag parameters whose values are given in the appendix. By setting C 
instead, one finds similar relations between f3 i and C . It follows that numerically 



(~) 
Pi 



>r>— ~tPf 
where we have used B. 



1.10 x 1(T 3 C g GeV 2 



P = -3.49 x 10 



-3 



C g GeV^ 



(17) 
= 

(18) 



m K/( m s + ™*)> witn m K = 495.7MeV and [2] m s + m = 121 MeV. 
We note that C and C here are the Wilson coefficients at the low scale /i — 0(1 GeV) and 



o 
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hence already contain the QCD running from the new-physics scales. We also note that the bag- 
model numbers in Eq. (18) are comparable in magnitude to the natural values of the parameters as 
obtained from naive dimensional analysis [26], e.g., 

flP 0024C » GeV2 • °- 0014C, « Gey2 • < 19 > 

where A = Anf is the chiral- symmetry breaking scale, with / = f n ~ 92.4 MeV, and we have 
chosen g s = \^4tc. The differences between the two sets of numbers provide an indication of the 
level of uncertainty in estimating the matrix elements. This will be taken into account in our results 
below. 

Next, we adopt the usual prescription for obtaining a weak phase [8, 12, 15], namely dividing 
the imaginary part of the theoretical amplitude by the real part of the amplitude extracted from 
experiment under the assumption of no CP violation. For the real part of amplitudes, we employ 
the experimental values obtained in Ref. [27]. They are, in units of G F m^. + , 

Sa-^t- = !- 42 ± 0.01 , s 3 .^ An . = -1.98 ± 0.01 , 

(20) 

p A _ p7T - = 0.52 ± 0.01 , P s -^ A7r - = 0.48 ± 0.02 . 

The imaginary part of the amplitudes are calculated from Eqs. (13) and (15), combined with 
Eq. (18). The other hadron masses that we employ are m N = 938.9, m A = 1115.7, m E = 1193.2, 
m s = 1318.1, and = 137.3, all in units of MeV. In Table I, we have collected the results, divided 
by the central values of Eq. (20), in terms of C g and C g . We find that in each of the amplitudes the 
(3p terms are numerically larger than the (3^ F terms, by up to a factor of four, and both contribute 
with the same sign. 



TABLE I: Ratios of theoretical amplitude arising from chromomagnetic operators to experimental ampli- 
tude, for 5- and P-wave transitions. 



Decay mode 


ImS (GeV) 




lmp (GeV) 




s expt 




Pexpt 


A — > pir~ 


-2.2 x 10 5 Im(C fl - 




-2.6 x 10 5 Im(C7 9 + C g ) 


S~ — » A7T~ 


-1.9 x 10 5 Im(C fl - 




+1.1 x 10 5 lm(C g + C g ) 



Since the ratios in Table I follow from the leading-order amplitudes in x?T, the uncertainty of 
our prediction will come partly from our lack of knowledge about the higher-order terms which are 
presently incalculable. Various studies of hyperon processes in the context of xPT show that the 
leading nonanalytic contributions to amplitudes can be comparable to the lowest-order terms [21, 
22, 24, 27, 28]. We expect that a similar situation occurs here, and consequently we also expect the 
uncertainty due to the higher-order contributions to be comparable to our leading-order estimate. 
To reflect this, as well as the uncertainty in estimating the matrix elements above, we assign an 
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error of 200% to each of these ratios. In Table II, we have listed the ratios as the weak phases, 2 
along with their uncertainties, in terms of = C g ± C . Accordingly, C g and C~ correspond to 
parity-even and parity-odd transitions, respectively. 

TABLE II: Weak phases generated by chrome-magnetic operators and strong-phase differences, 5 S — S p . 



Decay mode 


10" 5 4> s (GeV) 


10" 5 (j) P (GeV) 


S s - S p 


A — ► pir~ 


(-2.2 ±4.4) ImC; 


(-2.6±5.2)ImC+ 


7° ±2° 


S~ — > Air~ 


(-1.9 ±3.8) ImC; 


(l.l±2.2)ImC+ 


-2° ± 6° 



In Table II, we have also included the strong-phase differences. The number for A — > pn~ 
results from the measured phases quoted in Sec. I. For H~ — ► An~ , while awaiting a definitive 
measurement by HyperCP, we have adopted the range —7.8° < 5$ — 5p < +3.9° estimated in 
Ref. [17]. This range is compatible with the experimental values known to date, including the 
preliminary measurement by HyperCP mentioned in Sec. I. 

From the results in Table II, it follows that the contributions of the CMO are 
1(T 4 (A A ) = (3.5 ± 7.0) lmC~ + (-4.2 ± 8.3) ImC+ , 

(21) 

!0~ 4 (A s ) g = (-2.0 ± 6.0) ImC" + (-1.2 ± 3.4) ImC+ . 

where the numbers on the right-hand sides are all in units of GeV. This indicates that (A E ) g is not 
negligible compared to (A A ) and hence should be included in evaluating A EA . Summing the two 
asymmetries then yields 

ID" 4 (A SA ) g = (2 ± 13) ImC; + (-5 ± 12) ImC/ , (22) 

the right-hand side being again in GeV. The errors we quote here are obviously not Gaussian, and 
simply indicate the ranges resulting from our calculation of the phases. 

Since the CMO also contribute to the parameters e 1 and e in the kaon sector, it is possible to 
obtain a bound on their contribution to A EA using the measured values of e' and e. The contribution 
to e' can be written as [15, 18] 

= ^5.2 x 10 5 GeV) B G ImC; , (23) 

9 



2 We remark here that the central values of the numerical factors in front of Im for <pg P in Table II are larger, 
by roughly a factor of two, than the corresponding numbers obtained in Ref. [15], which considers contributions 
from a generic supersymmetric model. The disagreement may be due mainly to a factor-of-two difference between 
the matrix elements in Eqs. (Bl) and (B2) and those employed in Ref. [15]. 
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where B G parameterizes the hadronic uncertainty, and m s + rh = 121 MeV has been used [2]. The 
contribution to e occurs through long-distance effects, and the simplest ones arise from n , r], and 
rf poles [25], yielding 

(e) g = -(2.3 x 10 5 GeV) « ImC+ , (24) 

where k quantifies the contributions of the different poles. Hence Eq. (22) can be rewritten as 

. 0.04 ± 0.25 fe'\ 0.22 ± 52 . . 



To estimate the range of (A SA ) g allowed by the experimental values |e| = (22.80 ± 0.13) x 10 4 
and Re(e'/e) = (16.6 ± 1.6) x 10" 4 [2, 7], we require 



< 19 x 10~ 4 , |eL < 23 x 10" 4 . (26) 



Consequently, adopting 0.5 < B G < 2 and 0.2 < \k\ < 1, after Ref. [15], we find the bound 

|A SA | 9 < 97 x 10~ 4 . (27) 



The upper limit of this range is allowed by the published data [7] , but is disfavored by the preliminary 
result of HyperCP [6] quoted in Sec. I, exceeding it by several sigmas. Since the number in Eq. (27) 
is dominated by the (e) g bound, we can then conclude that the available preliminary measurement 
by HyperCP already probes the parity-even contributions better than e does. 

Now, it is possible that HyperCP will in the end observe no CP-violation in A 3A . In that event, 
the data can be used to estimate the bounds on both ImC^. To explore this possibility, we assume 
that HyperCP will be able to reach the expected sensitivity of 2 x 10~ 4 [6], and so we take this 
number as the upper limit for A BA . Moreover, since our result in Eq. (22) has large uncertainties, 
for illustrative purposes we use its central value in what follows. Barring significant cancellations 
between the ImC^ terms, we can consider three possible cases, (i) ImC^ = and ImC~ 7^ 0, 
(ii) ImC+ 7^ and ImC~ = 0, and (iii) ImC+ ~ — ImC~ 7^ 0. Consequently, requiring |A 3A | < 
2 x 10~ 4 , we obtain for these cases 



(28) 



(i) |lmC-| < 1 x lO^GeV" 1 , (ii) |lmC+| <4x lO^GeV" 1 , 

(iii) |lmC+| ~ |lmC-| < 3 x lO^GeV" 1 . 

For comparison, the requirements in Eq. (26) from e' and e measurements imply 

|lmC-| < 7.4 x lO^GeV- 1 , \ lm Cg\ < 5-0 x lO^GeV" 1 . (29) 

Similar or stronger bounds on ImC^ may also be obtainable from K — > 3ir, n£ + £~, 7177 decays [20]. 
Therefore, even if it turns out that HyperCP eventually does not discover CP violation in hyperon 
decays, its data can be expected to provide stringent constraints on the coefficients of the CMO in 
new-physics models, at a level that is comparable to or better than the bounds coming from the 
kaon sector. 
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IV. CONCLUSION 



We have evaluated the CP-violating asymmetries A A and A E induced by the chromomagnetic- 
penguin operators, whose coefficients can be enhanced in some scenarios of new physics. Including 
recent information on the strong phases in S — > An and adopting a model-independent approach, 
we have shown that (A s ) g , which was neglected in earlier studies, can be comparable to (A A ) g . 
We have found that the upper limit of the sum of these asymmetries, (A SA j g , as allowed by e and 
e' data is already disfavored by the preliminary measurement of A SA by HyperCP. It follows that 
this preliminary data already imposes a constraint on the parity-even contributions of the operators 
that is stronger than the bound obtained from e in kaon mixing. We have explored how well the 
upcoming results from HyperCP may bound the coefficients of the operators in the event of null 
results. In that case, the data will likely yield significant constraints that are comparable to or 
better than those provided by kaon measurements. 
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APPENDIX A: CHIRAL REALIZATION OF (3 L ,3 R ) AND (3 L ,3 R ) OPERATORS 

The form of the weak Lagrangian in Eq. (9) can be inferred from the lowest-order chiral realization 
of the operators d(l + 7 5 )s and d(l — 7 5 )s. The reason is that these densities, like the operators 
Q and Q g in Eq. (5), transform as (3 L ,3 R ) and (3 L ,3 R ), respectively, under SU(3) L xSU(3) R 

rotations. Using Cf^ in Eq. (8), one can derive the correspondences [12] 

-d L s R <=► b D {t}B v B v t} + ^B V B V ^) 32 + b F {$B V B V ? - ?B v B v t% + b Sj 2 (B V B V ) 

+ if 2 BA, (Al) 

-d R s L <=► b D (iB v B v i + iB v B v i) z2 + b F (tB v B v i-iB v B v £) Z2 + b Z 32 (B v B v ) 

+ \f 2 B ^ 2l (A2) 

where q L — |(1 — 7 5 )g and q R = |(1 + 7 5 )g. The form in Eq. (9) then follows. 

It is worth noting here that each of the S- and P-wave amplitudes in Eqs. (13) and (15) vanishes 
if we set 

c ~ ~ c 

Pd,F = C ^D,F 5 Pip = 2 ' Pd,F = C ^D,F 1 Pip = 2 ' (A3) 

with c and c being constants, and also use the relations in Eq. (14) as well as 

m E — m N = 2(b D — bp){m s — m) , m 5 — m s = —2(b D + bp){m s — m) , (A4) 
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derived from Eq. (8). This satisfies the requirement implied by the Feinberg-Kabir- Weinberg the- 
orem [29] that the operators d(l ± 7 5 )s cannot contribute to physical amplitudes [30], and thus 
serves as a check for the formulas in Eqs. (13) and (15). 



APPENDIX B: BAG-MODEL PARAMETERS 



Here we provide the estimate in the MIT bag model of the matrix elements of the chromomagnetic 
operators contained in Eq. (16). The relevant calculations can be found in Refs. [31, 32]. We have 

(n\g s da^(l± l5 )sG^\A) = ~ 16V *f» N4lM u n u A , 

(Bl) 

<A|,>-A«(1± 75 ) S ^|~°) = M 3 , 



(n-\g s da^X a (l± l5 )sG^\K-) = ^M^M (B2 ) 

where A a = 2t a , only the parity-conserving part of the s — > dg operators contributes, and R, 
N, and I M are bag-model parameters [32]. Numerically, we choose g% = 4n, corresponding to 
a s = 1, and adopt R = 5.0 GeV -1 for the baryons and R = 3.3 GeV -1 for the mesons [32]. 
Since the weak parameters f3 i and f3 i belong to a Lagrangian which respects SU(3) symmetry 
[C w in Eq. (9)], in writing Eqs. (Bl) and (B2) we have employed SU(3)-symmetric kinematics. 
Accordingly, we take m u = m d = m s = and use the formulas given in Ref. [32] to obtain 
N = 2.27 and I M = 1.63 x 10~ 3 for both the baryons and mesons. Finally, we note that Eq. (B2), 
together with the relation (ir~\Q g \K~) = — {ir°\Q g \K°), leads to the matrix element [25, 33] 
A Rn = <7r°|^A a (l + 75 )sG^ ) = +64 g s N 4 I M m K /R 2 ~ 0.4GeV 3 . 
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